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ABSTRACT 


A truncated procedure for selecting the better of 
two binomial populations is formulated. The procedure 
uses play the winner rule and the selection is based on 
the difference in the proportion of successes at termina- 
tion. The trials are terminated whenever either population 
is sampled completely or when the difference in the 
proportion of successes exceeds a critical value, based 
Onethe taliures on both populations. 10 1s shown that. 
when the average of the two population probabilities of 
Success is less than 0.5, the expected number of observa- 
tions on the poorer treatment is smaller for this proce- 


dure than for any procedure currently known. 
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CHAPTER I 


INTRODUCTION 


The theory of clintcal trials is of considerable 
interest to theoretical and applied statisticians. The 
idea behind most of the newer techniques is the ethical 
One of not prolonging a trial longer than necessary, for 
a trial which is unduly prolonged may result in an exces- 
Sive number of patients being given the less beneficial 
treatment. The doctor treating a patient in a clinical 
trial is not only obliged to derive information about 
the best treatment, but is also obliged to treat each 
patient in the best way possible. This ethical considera- 
tion has motivated the development of several statistical 
techniques for selecting the best treatment, or equiva- 
tently the selectton of the best of several binomials. 

The past work in selecting the better of two 
binomial populations falls into three categories: 

(a) aspects of the allocation problem, (b) two-armed- 
bandit problem and (c) alternate ways of conducting 
clinical trials by adopting the sampling technique. 

In the two-armed-bandit problem, we deal with a 


slot machine. When the left arm is pulled, the machine 
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pays one unit with probability p, while a pull on 

the right arm pays one unit with probability p'. 

The problem with a two-armed bandit is which arm to 
pull. Assuming the probabilities of pay-off of the 
two arms are constant, a strategy which will maximize 
the expected winnings for a sequence of pulls*is 
desired. 

In the allocation problem, the total number of 
patients waiting for the treatment is known. Each 
patient is to be treated with one of the two treatments. 
The first n patients are used to select the better 
drug and the remaining patients are treated with the one 
selected as being better. 

In practice, no knowledge about the number of 
patients is available. Zelen[11] modified the allocation 
model by assuming that (a) two (or more) similar groups 
are to be observed where the groups differ only in the 
treatments received; (b) the construction of the groups 
is by random allocation; (c) the response is dichotomous 
and the results are observed without delay. He suggested 
the use of sampling techniques to reduce the number of 
patients put on the poorer treatment. Zelen applied 
Play-the-Winner(PW) rule to the clinical trials; it 
prescribes that a success with a given treatment generates 
a further trial on the same treatment while a failure 


generates a trial on the alternative treatment.This is 
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a very simple rule to use that introduces bias in 
favor of testing the better treatment. 

Sobel and Weiss[9] compared different procedures 
for selecting the better treatment. A procedure is 
specified by the sampling and termination rule. A 
comparison of PW rule was made with Vector-at-a-time(VT) 
rures Vl consists Of taking two observatrons at each 
stage, one from each population and one does not consider 
stopping between these two trials. The situation where 
a physician wishes to compare two or more treatments 
using the patients he ordinarily encounters, and the 
problem of selecting the most effective dose of a 
particular drug are two among the several practical 
situations where these simple methods can be applied. 
Most of the investigations are done by formulating the 
selection of the better of two binomial populations 
(i.e. the one with the highest probability of success 
p on a single trial) as follows: for preassigned constants 
p’ and ,* with 0.5 Daten Mand an) a4 enc 
required that the probability of correct selection(CS) 
be at least P*when the true difference in p-values 
(denoted by, ) is at leastaAx . 1.€. we want a procedure 
R such that 

PCCS7 Rie p™ whenever he ae ahs 
Among the procedures satisfying this requirement, the one 


with the smallest expected number of patients on the less 
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beneficial treatment is declared to be the better 
treatment. 

In chapter II, a review of the recent work done 
on the binomial selection problem is made. In chapter III, 
we discuss a procedure R(DPPW) using PW sampling where the 
rule is based on the difference in the proportion 
of successes. The performance of this procedure is 


evaluated in chapter IV. 
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CHAPTER 1 


A REVIEW OF THE RECENT RESULTS ON 
BINOMIAL SELECTION PROBLEM 


In this chapter, we will briefly review the recent 
work on selecting the better of two binomial populations. 
The results mentioned in this chapter will be used for 
comparison with the results of the next chapter. Let us 


denote by 


S(i), F(i), and N(i) the number of successes, failures 


and trials respectively on the ath 


population for i = A,B 
N = N(A) + N(B), the total number of observations on 
both populations. 
A = p - p' where p and p' are the probabilities of 
success on a single trial on population A and B respec- 


Liveryer 16t, 0 25) 
CS - correct selection, namely identifying population 
A as the better one. 
Coy 
( P , aA* ) - preassigned constants such that 


Grew P xa and (Pee? oO 
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We need a procedure R for selecting the population with 


the larger probability of success on a single trial such 


that 


Goo PL CS/ Rt ieee whenever 4 > a*., 

Sobel and Weiss [9] considered procedures R[DPW ] 
and RLDVT] with the termination rule based on the difference 
in successes at each stage. The procedure R[DPW] uses play- 
the-winner ( PW ) sampling and declares treatment i as the 
better one when S(i) - S(j) = r, where j is the other treat- 
ment, while R[DVT] uses vector-at-a-time ( VT ) sampling 
and makes the same decision when S(i) - S(j) = s, The 
integers r and s are chosen to be the smallest such that 
(2.2) is satisfied for preassigned constants fone A* ) 

Sac ust ying (2st). 
The probability of correct selection under R[DPW] 


is calculated by defining 


Pine = cP GS / Si Ajo SB) =ene. Nt ans) 
C28 ) 


O(n) = PCS SCA) 1 seGb Ria T= Bo 
where NT = A (NT = B) denotes the next treatment is on A(B). 
Since the first treatment is chosen at random, 


the P(CS) at termination is given by 


— 


(2.4) P(CS/RLDPW]). <p py Pi Ore D0) J 


Onder RI DPW), letting g = 1 o= ps andq = l= p°s Wwesobtain 
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with the boundary conditions» P(r) =*1.°0(=r) = 0% 

The solution to (2.5) satisfying the boundary conditions 
is given by 

rH ee q'=gig " 

where A= (p'/p)< 1. Setting P(CS)=P., we obtain from (2.6) 
ange cc. 4) 


2-160) Q(n) = 


C2 q'- 5 (qtq')a® = P(q'-qa2") 
as the equation determining r. Thus given any p and p', 
(2.7) determines r such that PYCS) =P. Since it is desired 
that (2.2) be satisfied when p and p' are unknown, it 
is necessary to find the Least Favourable ( LF) configura- 
tion ofp and p* which Sadtisfies (2.2). By LF we mean p 
and p' which . minimizes P(CS) subject to the restriction 
p-p'> A*. It is shown in [9] that, when cae 1, the value 
of r under LF can be approximated by 
Ee 
r = ee eee een 
" 1og(1-a*) 
where [.] denotes 'the largest integer contained in’. 
Under the procedure R[DVT], using VT sampling, we 
stop when S(i)-S(j)=s and choose the treatment i as the 


better. We let Pin) = P(CS/S(A)=S(8) =n)e ) thus 
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P(CS/R[DVT]) = P(0); and it is obtained by solving 


the difference equation 
( 2143)) P(n) = pq'P(nt+1)+qp'P(n-1)+(pp'taq')P(n) 


With the boundary conditions P(s)= 1 and Pi-s)=0 
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It is shown that minP(CS/R[ DVT) for A = A* is attained 
bywsetthing Ay = A* “and? p= (A+) /2.. The requireds 


is the solution of 


(2.9) 


* 
It is shown in [9] that, when P + 1, fhe value of s under - 


LF can be approximated bg P 
Slogtl=P. 


The procedures R[DPW] and R[DVT] are compared by the expect- 


ed number of trials on the poorer treatment. Defining 
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then E(N(B)/R[DPW]) = 5 R(0)+S(0)} 
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Under R[DPW] we obtain 


iT] 


RGn) = pR( ntl rary gS.(n.) 


esti.) 
Se = pS ina eq k(n) 4 


with the boundary conditions R(r) = S(-r) = 0 


Solving (25512). We obtain 


+2 Ten (6 oe Ie 
E(N(B)/ ppw) = L (ee | 
q'-qa* 


where A = p'/p . Under R[ DVT], we define 


(21.3) G(n) = E(N(B)/ S(A)-S(B) = n ) 


we obtain therefore 


(2.14) a(n pao Gint)) Gp Gi no pp ageaia cht 
with the boundary conditions G(s) * G(-s) = 0. Solving 
(2.14) we obtain 


-§(1-8°) 
E(N(B)/R[DVT]) = + {| 


where 6 = p'q/ pq' . Further it was shown in [9] that for 

on close to 1 and A* small, R[ DPW] is preferable to 

R[DVT] when p> , A - £-} 5 otherwise R[DVT] js 

preferred. Numerical results of E(N(B)) for different 

(P- a*) combinations are given in Tables 6.) 75.8 and: 9); 
Sobel and Weiss [8] considered an inverse stop- 

ping rule which terminates the. trials.when any one popu- 


lation has r successes. The procedure R[ISPW] uses 
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PW sampling with inverse Sampling as the termination 


rule. It was shown in [8] that 


P(CS/R[ISPW]) = 5 J ctin ail Gielen 


where 


and ELf(x),| is the expectation of f(X%) where Xx 1s a 
(negative binomial) random variable denoting the number 


th 


of failures before the r success. A normal approximation 


to minimum of P(CS/R[ISPW]) when r > » is given by 


(2.15) Min P(CS/RLISPW]) ey aly for AzA* 
8 


where (x) is the standard normal c.d.f. We solve for r by 
putting A=A* in (2.15) and setting the result equal to 
p>. if oe »(P”) denotes the solution of o(\) =P , 
then we obtain 
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An asymptotic ( r+ ~ ) normal approximation to E[N(B)] 


under R[LISPW] is given by 


E(N(B)/RCISPW]) = ©, ¢ 2 aly) + SL ( 1-0(y))) 
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where y =a | = , D = aq(p')® + q'(p)2 ° As seen in 
tables 6,758 and 9, the inverse sampling rule requires a 


large number of observations on the poorer treatment 
when the success probabilities are small. Another 
procedure RLIFPW] resulting in smaller E(N(B)) when the 
success probabilities are small, is obtained by changing 
the termination rule so that the experimenter will wait 
T Ota. Fixed number of failures instead of successes. 
Under R[IFPW], one waits for r failures from each popu- 
lation and selects the one with most successes as being 
best, using randomization when there is a tie. It was 
Shown that P(CS/RLIFPW) = P(CS/RLISPW]) and the critical 
values r are given by (2.16). It was further shown, when 
reo, E(N(B)/RLIFPW]) = r¥/q~ @ tables 6 74.0 ana. 
give E(N(B)) for different (P >, a*) combinations. 

Hoel [3] modified RLISPW] to obtain a truncated 


test, R[H], that uses PW sampling. Defining 
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R[H] terminates the trials as soon as R(i) or R(j) is 
equal to r, where r is chosen such Chat ceee el oS aS ted 
If R(i) equals r, then population i is declared to be the 


better one. The LF configuration was found to be 
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P(CS/R[H]).. = Eveq t Ti (Xt, r-X) ] 


where ace) denotes the expectation of the binomial 
random variable X with parameters r and q and I, (b,c) es 
the incomplete beta function. We find E(N(B)/R[H]) by 


defining 
pee BNR OAD tM re RD ae 


RCA RUBY ase im= (2.16) 


R(m,n) ELN(B)/ T=(m,n), NT=A] 


S(m,n) = ELN(B)/ T=(m,n), NT=B] 


Under R{H] we obtain 


R(m,n) p R(m-1,n) + q S(m,n-1) 
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Tables 6, 7, 8 and 9 show that R[H] has a larger expected 


sample size on population B than RL[LISPW] whenever 
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Another truncated procedure in which the maximum 
number of tests was specified was discussed by Kiefer 
and Weiss [6]. This procedure uses VT sampling. After 
each pair of populations is sampled, one calculates 
MSV S(i)=SCj ce Th at any test ti< NO the maximum 
number of tests) we have AS = s, then we terminate, 
calling the population with the greater number of suc- 
cesses the better population. If no decision is made 
within the N test pairs, the populations are held to 
be essentially of equal value. In this truncated version 
the experimenter can use the parameter N as follows: 

(a) Fix N and determine s so that P(CS) is maximum. 

(b) Choose the smallest N consistent with (2.2). 

(c) Require, in addition to (2.2), that the proba- 


* 
bility of a decision of equality be = P, when 


The authors were not successful in deriving an exact 
expression for the expected number of trials to reach a 
decision. Values*of s“that maximize P(CS)*tor fixed a* 
and N are given in [6]. 

Hoel, Sobel and Weiss [4] discussed a two-stage 
Sampling procedure: In thegfirst stage \N tests are cones 


ducted on each population. The type of sampling to be 


used in the second stage was based on a critical parameter 


k ; PW sampling was used in the second stage if the 
maximum of the number of successes in the first stage 


was > k ; otherwise vector-at-a-time sampling (VT) 
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was used. The terminal decision rule was based on the 
difference in the successes on both populations. Expected 
sample sizes on the poorer population for different (p>, a*) 
combinations are given in chapter IV. It was shown that 
the two-stage procedure resulted in substantial savings 
On the expected samples on the poorer population compared 
to R[DVT] whenever p > 0.6. We denote this procedure by 
REDTS.) 

Nebenzahl and Sobel [7] investigated the selection 
problem for fixed sample size case. The procedure R[FSPW] 
uses PW sampling while R[FSVT] uses VT sampling. R[FSVT] 


forces the sample size to be an even number. The same 


Table 1. 
Comparison of R[FSPW] and R[FSVT] for the pair 
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terminal decision rule is used for both procedures ,namely 
to select the population with most successes and to rando- 
mise when we get equality in the number of successes. The 
fixed sample size N for both procedures is determined so 
that (2.2) is satisfied. Table 1 gives N needed under 
REFSPW] and E[LN(B)] under R[EFSPW] and REFSVT] for differ- 
ent ee A*) combinations. 

Another procedure R[IT] due to Berry and Sobel [2] 
modifies the inverse sampling procedure R[ISPW] by termi- 
nating the trials either after c failures On each population 
or when there are r successes on any population, whichever 
occurs sooner. In either case the population with the 
larger number of successes is selected. The E[N(B)] when 
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CHAPTER ITI 


A TRUNCATED PROCEDURE FOR CHOOSING THE 
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The different procedures discussed earlier are 
distinguished by the sampling and termination rules. 
The termination rule was specified by a decision 
function, defined by the current number of successes 
andor twan lures on either or both populations, and a 
critical value r, a preassigned constant chosen to 
Satisfy (2.2). In this chapter we analyse a truncated 
procedure using PW sampling, in which a maximum number 
of tests on each population is specified. 

Consider two bernoulli populations, one having 
Orobabillity Of Success D on acsingle trial and probability 
of failure q =1]- p; denote it by A; the other is B 
with parameters p'(p>p') and g'=]-p' .The parameters 
p and p' are both unknown and our sampling and termination 
rule will not depend on the knowledge of which population 
isehaent the oubuset one Of the two populations 1s enosen 
at random. Procedure R[DPPW] uses the PW sampling with 
the selection based on the difference in the proportion 
Of Successes at termination. Before the trials are 


started, the maximum number of tests, N, on each population 


Slee 


ore Aekithea koedusebs eaauhaneeg dnaret NaN wir) oat / 

,2oToy fokisenimist bas prt lqamsa ons XG | 

not2toeb & Vd bet ttoaqe 25w 3fiy notsentmrssoedT 

eseeedsue to Vedmuq srgyswo ont yd Bont ted snot soma? a) 

b bos. -enofselugeg Atod yo VAT hse sowie? vo\bas 7 

ot nezada Intends benptzzserq 6.1 ovley Ieottiia 

bstsonuad « seylang aw totaads ett wl «{8.8) qhettee | 

jedmun mumtxsm 6 dotdw at ,entiqmee WS gniew oF nq 7 : 

bat ttoege 2f mortsfugoq dase ao 2deed 20 a 

Prive ano g2netsniygqaa if ueored ows Ysbtened Z 

ysl tdedong bas Taiwt elpate 6 mo 9 e2so0u2 to on 
a 2+ sedto sd GA Yd th stonsb gq -7* p owrte? to 

27S F3tbVEq gut. gq= P= 'p bne (iq <q) 'q 219t9m615q wt 

rortsAiiyes bre ent Tqmee 100 BE AWondny tod Svs "q bs q 

alt Rac tah att wo hesqsb ton Tfiw stun 

2h z2nohtufugog owt emt Yo sno dscdue add 2h WA et. a 

acute at caru [W99G]9 Syiban0"8 -mobaey $6, 

at | aid no baz notszeise ey | 

ana SYote® <notsentmiet 25 cozesague Lee 
aan vod anwn mumixem ott yt J 


an 
“at a 


2S 


ond c Lhevcritical parameter are specified. Using 

the notation in chapter II, the procedure R[DPPW] has 
the following stopping rule: sampling is terminated 
whenever either population has been sampled completely 


or when 


N(A) N(B) GAY FB) 
whichever is earlier, where the constants (N,c) are 
chosen to be the smallest integers satisfying the 


condition 


* 


(3.2) P(CS/R[DPPW]) > P whenever 4 > q 
for preassigned constants crops AX) Satrstytng. C2.2)e 
The population with the larger proportion of successes 


at termination is declared as the better. This procedure 


is truncated, since at most 2N-1 trials will be required. 


The results of the total sample size required for the 
fixed sample size binomial selection procedure R[FSPW] 
given in Table 2 were used as follows to specify (N,c). 
For any specified (PF » Ax) combination, N was taken to 
bé slightly greater than half the total sample in the 
fixed sample procedure with 

fees) P(CS/R[DPPW])> P when 4 = 4* 
Then, for this N, the smallest tmveder cesatist ying 


(20) was found? then Skor This yale yortet ithe smallest 
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N satisfying (3.3) was found. The combination (N,c) found 


thus was taken to be the final choice of (N,c). 


bape 2 


Total Samples N required under R[FSPW] 
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The expected sample sizes on the poorer population, when 
the population probability of successes are small, are 
considerably reduced by this procedure. 

To Galicuiater P(CS/REDPPWLI.. bet 

NT=A(NT=B) denote the next trial is on population A(B). 


T = (m,n,a,b) where 


m= N - N(A) 

n = N - N(B) 

a = number of failures on A 
an) 


BD Einumperk Ota failures on B 


Define 


Utmsn,a,b) = P(CS/, J=(m.n,4.5)  NT=A) 
(35) 
Vinemaasb) <= P(CS/ =(mnsaabjiaiNi=B) 


Using PW sampling and using [ (3.4),(3.5)] we obtain 
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for the procedure R[DPPW], 


U(m.n,a,b) = p U(m-1,n,a,b) + q V(m-1,n,at+1,b) 
(3.6) 


Vim,n,a,b) = p'V(m,n-1,a,b) + q'U{(m,n-1,a,b+1) 


and the boundary conditions are 
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(gy Oa ee eer es Dene Ne [a-b] < 1 


The first treatment at the’ outset. is’ chosen at 


random and hence P(CS) at termination is given by 


(3.8) P(cs/rtpppw]) = + (u(n,n,o,0)+v(N,N,0,0) | 

An analytical solution to (3.6) with the boundary 
conditions(3.7) is quite complicated. We used the IBM 360 
computer to numerically evaluate the P(CS/R[DPPW]) 
when N and c are specified. Several computations using 
program A (appendix) were made to find the value of c 
and LF configuration for preassigned (N, as and - 


the results are summarized below. 


Table 3. 


The Least Value of c for R[DPPW] satisfying 


(P's N) probability requirement for specified N 


The performance of R[DPPW] was evaluated by comparing 


the expected number of observations on the population B 
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Program B (appendix) was used to evaluate 
E(N(B)/R[DPPW]) for different qe, A*) combinations and 


the results are summarized below: 


Table 4 


Expected Number of observations under R[DPPW] 
ELN(B)] : Expected number of observations 


One POPuUl akon) iB. 


ELN)| : Expected total number of observations. 
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required for P =,.,90 and A =.1 


Consider the procedure R[MDPPW] obtained by 
modifying the stopping rule of R[DPPW]. The-procedure 
RLEMDPPW] uses PW rule and sampling is terminated 


whenever either population is sampled completely or 


s(A) $(B) | cf S(A)#8(B) 
N(A) — N(B) F(A)+F(B) l N | 


when 


where [ .] denotes " largest integer contained in". 


It was verified numerically when c=5, N=38, the 
P(CS/RLMDPPW]) to be at least .95 whenever p - p's.2. 
Wieelr contigurattonsis: given by (ip,p )) = (75208 32) 
and the expected sample size on the poorer population 
(B) is summarised in Table 5. The computations were 
made using program A and program B listed in the 


appendix. 


It should be noted that R[MDPPW] attains 0.95 


probability of correct selection with a smaller N and 
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Table 5. 


Comparison of E(N(B) for the procedures 


REDPPW] and R[MDPPW] for (P ,4) =(.95,.2) 


performs uniformly better than R[DPPW]. The procedure 
REMDPPW] has not been investigated thoroughly. In the 

next chapter, we discuss the unsolved problems | 
regarding R[DPPW]. The performance of REDPPW] with R[ISPW], 


RL IFPW],R[DPW],REDVT],RCH],R[DTS] and R[IT]Jdiscussed earlier. 
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CHAPTER IV 
COMPARISON OF R[DPPW] WITH OTHER PROCEDURES 


Procedure R is said to be uniformly better 
than procedure R' if E[N(B)/R] < E[N(B)/R'] . In 
this chapter we use this criterion to compare the 
different procedures described in the previous 
chapters. Later we state some open questions concem- 
ing the procedure R[DPPW] and some open problems 


in selecting better of two binomial populations. 


RLDPPW] is a truncated sampling procedure with 
the selection based on the dtfference in the proportion 
of successes at termination. Tables 6, 7, 8 and 9 
give the expected number of observations on the poorer 


population for various procedures. Let p = (ptp')/2 
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Table 6. 
Expected number of observations on the 
poorer population. 
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Table 7. 


Expected number of observations on the 


poorer population. 
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Table 8. 


Expected number of observations on the 
poorer population. 
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Table 9. 


Expected number of observations on the 
poorer population. 
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Table 10. 


Expected number of observations on the 


* * 
poorer population: SCP. Ay \eiarGeo 5s 2) 
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The following conclusions are made from the numerical 
Yeic Wile S) GiNus lidib |eSu0.7 2 oad Ge or 
(a) We observe that whenever p < 0.5, the procedure 
RLDPPW] is uniformly better than any other pro- 
cedure currently known; R[DPPW] results in a 
Substantial more than 25%) Savings in ithe 
expected number of observations on the poorer 
population. 
(b) We conclude that for p > 0.7, R[DPW] is uniformly 
better than any other procedure. 
(c) When 0.5 < p < 0.7, none of the procedures is 
uniformly better. However, if the experimenter 


— * 
knows apriori that p = 0.6, then for P = 0.90, 


A = A*= 0425 RIOTS) is recommended. 


Table 10 gives the expected number of observations 
on the poorer population when A > A*. The performance of 
R[DPPW] is compared with R[H], R[DPW] and R[DVT]. As 
expected for p < 0.5, the performance of R[DPPW] is 
better than any other procedure. 

The better performance of R[DPPW], when p its small, 
can be intutively justified for the following reasons. 
Since the maximum sample size, N, is specified in advance 
R[EDPPW] prevents unduly long continuation of trials. The 
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parameter,c, and the total number of failures at each 
stage. The critical value decreases whenever a failure 
OCCUNSIZ witch: ths; contrary eto sthhe<constan-t, .cimi.tica ly.v.alue 
specified by the other procedures. If p < 0.5, the result 
is frequent occurence of failures, and termination will 
occur when the difference in the sample proportion of 
Successes exceeds the gradually decreasing critical value 
when N is large. However if N is small, we terminate the 
trials by sampling either population completely. In either 
case, we expect R[DPPW] to terminate the trials quicker 
than other procedures, due to the decreasing characteristic 
of the critical value. Consequently, we expect R[DPPW] to 
discontinue sampling with a smaller total number of fail- 
ures. The procedure is specified such that fewer failures 
are needed to distinguish a large difference in the sample 
proportion of successes.For example when P~=0.95, c=5 and 
A Fabrice 0 2s, R[DPPW] requires only 13(26) failures to 
terminate the trials when the difference in the sample 
proportion of successes 45° 0.2(0:1). Further the critical 
parameter prevents us from terminating trials at a very 
early stage, thereby reducing the probability of making a 
wrong selection by chance at a very early stage. If p> 0.5 
the result is frequent occurence of successes and the criti- 


cal value does not decrease quickly and hence we terminate 


the trials by sampling either population completely. 
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Consequently, in this situation the performance of R[DPPW] 
is not better than some procedures. 


The principal shortcomings of the procedures described 
above are the assumptions of dichotomous and instantaneous 
response to treatment. So far, no analogue has been pro- 
posed for the continuous case in which the response is 
continuous - a very important consideration in testing 
anti-cancer treatments for which a natural measure of 
effectiveness is life time. Further this problem has not 
been analysed when the response is at some random time 
after the administration of the treatment. The selection 


problem under the present formulation has not been studied 


when several patients are assigned to any treatment at each 


Stage. In fact, any enumeration of practical difficulties 
associated with clinhedl trialist is tasgood source for 
future investigation. Most of the stopping rules excepting 


the inverse sampling lead to mathematical difficulties. 
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DIMENSION U (66,66,3) ,U1(66, 66,3) ,V (66,66,3) ,V1 (66,66, 3) 
THIS PROGRAM CALCULATES THE PROBABILITY OF CORRECT 
SELECTION. IN THIS PROGRAM 
MP 1= (M+1) 
NP 1= (N+1) 
_IFAP 1=A+1(#FAILURES ON POPULATION A + 1) 
IFBP1=B+tl (#FAILURES ON POPULATION B+ 1) 
IFAPFB= TOTAL NUMBER OF FAILURES 
WHERE (M,N,A,B) ARE DEFINED BY (3.4). P(P1) IS THE 
PROBABILITY OF SUCCESS IN A SINGLE TRIAL ON POPULATION 
A(B). DELy IS LHE DIRFERENGE IN -SUCCESS PROBABILITIES. 
IR IS THE MAXIMUM SAMPLE SIZE ON EACH POPULATION, 
THE PROBABILITY OF CORRECT SELECTION IS GIVEN BY (3.8). 
FOR THIS WE CALCULATE U(N,N,0,0) AND V(N,N,0,0). 
FIRST WE FIX MP1; AND THEN ALLOW NP1 TO VARY. NOW FOR 
FIXED NP1, WE LET IFAP1 TO VARY; SINCE PW SAMPLING IS 
_IS USED IFAP1 AND IFBP1 DIFFER BY AT HOST 1. IFBP1 IS 
DEFINED IN TERMS OF J3 WHERE J3=IPBP1-IFAP1+2. 
IN THIS PROGRAM EACH U(M,N,A,B) AND V(M,N,A,8) IS 
ASSOCIATED WITH A PARTICULAR VALUE OF MP1 AND 
U(NP1,IFAP1,J3) AND V(NP1,IFAP1,J3) RESPECTIVELY; 
BY FIXING MP1 FIRST, WE NEED TO USE ONLY THREE DIMENSIONS 
IN THE PROGRAM, THEREBY MAKING A SUBSTANTIAL SAVING OF 
VIRTUAL MEMORY STORAGE. THE VALUE OF U (M-1,N,A,B) 
NEEDED FOR CALCUATING U(M,N,A,B) IS STORED BY SETTING 
O1(NP1,IFAP1,J3) = U(NP1,I1FAP1,J3). 
WE FIRST LET U(M,N,A,B)=0, V({(M,N,A,8)=03 LATER WE 
CHECK THE BOUNDARY CONDITIONS AND TERMINATION 
CRITERION. : 
IF REQUIRED, U(M,N,A,B) AND V(M,N,A,B) ARE CALCUALTED 
BY USING THE DIFFERENCE EQUATIONS. 
READ(5,14) P,DEL,IR 
FORMAT (F10.6/F10.6/12) 
Q=1.-P 
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Q1=1.-P1 
IRP1=IR+1 
DO 15 MP1=1,IRP1 
DO 15 NP1=1,IRP1 
NAP 1=IR+2.-MP1 
DO 15 IFAP1=1,NAP1 
NBP1=IR+2-NP1 
WHEN J3=1, A=Btl; WHEN J3=2, A=B; WHEN J3=3, B=A+1. 
DO. 15.J03=1, 3 zo ; 
IFBP1=IFAP1+J3-2 
U(NP1, IFAP1,J3) =0 
V(NP1,IFAP1,J3)=0 
WE CHECK THE BOUNDARY CONDITIONS AND TERMINATION 
CRITERION 
IF(MP1.NE.1) GO TO 44 
MP1=NP1=1 MEANS THAT BOTH THE POPULATIONS ARE SIMULTANEOUSLY 
SAMPLED COMPLETELY, WHICH IS NOT POSSIBLE IN PW SAMPLING.. 
IP (NP1.E0.1)GO TO 16 
NUMBER OF FAILURES ON POPULATION B CANNOT EXCEED THE NUMBER 
TRIALS ON THAT POPULATION. 
IF(IFBP1.GT.NBP1) GO TO 16 
J3=1 MEANS A=B+1; THE NUMBER OF FAILURES ON POPULATION B 
IS ONE LESS THAN THE NUMBER OF FAILURES ON POPULATION A. 
SINCE THE NUMBER OF FAILURES IN PW SAMPLING DIFFER BY 
AT MOST 1, THE NEXT TRIAL CANNOT BE ON POPULATION A. 
IF (J3.EQ.1) GO TO 18 
IF(NP1.EQ.IRP1) GO TO 20 
IP(IFBP1.£Q.1) GO TO 18 
U(NP1,IFAP1,J3) =1 
IF (IFAP1.EQ.1) GO TO 16 
WHEN J3=3, B=A+1; SINCE THE NUMBER OF FAILURES IN PW SAMPLING 
DIFFER BY AT MOST 1, THE NEXT TRIAL CANNOT BE ON POPULATION B. 
IF (J3.EQ.3) GO TO 16 
IF (NP1.EQ.IRP1) GO TO 21 
IF( IFBP1.EQ.1) GO TO 16 
V(NP1,IFAP 1,33) =1 
GO TO 16 
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IF (IFBP1.GT.NBP1) GO TO 16 
IF(NP1.£Q.1) GO TO 16 
IFAPFB DENOTES THE TOTAL NUMBER OF FAILURES. 
IFAPFB=IFAP1+IFBP1-2. 
IF (IFAPFB.LE.1) GO TO 24 
RULE IS BASED ON THE TOTAL NUMBER OF FAILURES AND THE CRITICAL 
PARAMETER C; IT SPECIFIES THE CUTOFF POINT FOR THE TERMINATION 
CRITERION. 
RULE=5./IFAPFB 
DIFPRO IS THE DIFFERENCE IN THE PROPORTION OF SUCCESSES. 
DIF PRO=((IFBP1-1.) /(NBP1-1) )-—( (IFAP1-1.) / (NAP 1-1) ) 
DIFNEG=-DIFPRO 
IF(J3.EQ.1) GO TO 23 
IP(NP1.EQ.IRP1) GO TO 26 
IF (IFAPFB.LE.1) GO TO 26 
‘WE TERMINATE AND MAKE THE CORRECT SELECTION IPF 
THE DIFFERENCE IN THE PROPORTION OF SUCCESSES 
IS GREATER THAN OR EQUAL TO THE CRITICAL VALUE 
(RULE); OTHERWISE P(CS/T=(M,N,A,B) IS CALCULATED 
BY USING THE DIFFERENCE EQUATIONS. 
IF (DIFPRO. LT. RULE) GO TO 31 
U(NEP1,IFAP 1,33) =1 
GO TO 23 
WE TERMINATE AND MAKE A CORRECT SELECTION IF -DIFNEG IS LESS 
IS LESS THAN OR EQUAL TO RULE. 
IF(DIFNEG.GT.RULE) GO TO 23 
CALCULATING P(CS/T=(M,N,A,B) NT=A) BY USING THE DIPFERENCE EQUATION. 
U(NP1,1FAP1,J3) =P*U1(NP1,IFAP1,J33) +Q*V1(NP1,IFAP1+1,J3-1) cee 
IF (U(NP1,IFAP1,J33) .LT.10E-10) U(NP1,IFAP1,J3) =0 
IF(J3.EQ.3) GO TO 16 
IF(MP1.EQ.IRP1) GO TO 27 
IF (IFAP1.£0.1) GO TO 16 
IF(IFAPFB.LE.1) GO TO 27 
IF (DIFPRO.LT.RULE) GO TO 33 
V (NP 1, IFAP 1,33) =1 
GO TO 16 
CHECKING THE DIFFERENCE OF PROPORTION CONDITION. 
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IF(DIFNEG.GT.RULE) GO TO 16 


CALCULATING P(CS/T=(M,N,A,B) NT=B) BY USING THE DIFFERENCE EQUATION. 


V(NP1,IFAP 1,33) =P1*V1(NP1-1,I1PAP1,J3) +Q1*01(NP1-1,IFAP1,J3+1) 
IF(V(NP1,IFAP1,J33).LT.10E-10) V(NP1,IFAP1,J3) =0 
VALUE OF U(M-1,N-1,A,B) IS STORED FOR CALCULATING 

U(M,N,A,B) BY SETTING U1(NP1,IFAP1,J3)=U(NP1,IFAP1,J3) 


_U1(NP1,IFAP1,J3) =U (NP1, IFAP1,J3) 


V1(NP1,IPAP1,J3) =V (NP1,IFAP1,J3) 
CONTINUE 

CALCULATION OF PROBABILITY OF CORRECT SELECTION. 
PCS= (U1(IRP1,1,2) +V1(IRP1,1,2))/2 

WRITE(5,210) IR,U1(IRP1,1,2) ,V1(IRP1,1,2) ,PCS 
FORMAT (* *,*MAXIMUM SAMPLE SIZE',1I5,2X,'U =',F12.6, 
G2X,* V=',F12.6, "PROBABILITY CORRECT SELECTION',F12. 6) 
CONTINUE 
STOP 

END 
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DIMENSION R (73,73, 3) ¢R1(73,7343) 6S (73 673 63) o81(73473,3) 
THIS PROGRAM CALCULATES THE EXPECTED NUMBER OF 

PATIENTS ON THE POORER TREATMENT. IR DENOTES THE 

MAXIMUM SAMPLE SIZE AND DEL DENOTES THE DIFFERENCE 

IN SUCCESS PROBABILITIES. 

THE NOTATIONS AND STRUCTURE OF THIS PROGRAM ARE IDENTICAL 
TO PROGRAM A, THE BOUNDARY CONDITIONS ARE GIVEN BY (27s 
READ (5,14) P,DEL,IR 
FORMAT (F10.6/F10.6/12) 

FORMAT ('0°,T45,'MAXIMUM SAMPLE SIZE", 170,23) 

PRINT522,1R 

Q=1.-P 

P 1=P-DEL 

Q1=1.-P1 

IRP1=IR+1 

pO 15 MP1i=1,IRP1 

pO 15 NP1=1,IRP1 

NAP1=I1R+2.-MP1 

DO 15 IFAP1=1,NAP1 

NBP1=IR+2-NP1 

po 15 J3=1,3 

IFBP1=IFAP1+33-2 

R (NP 1, IFAP 1,33) =0. i 

S (NP1,IFAP1,J3) =0. 

IF((IFAP1.EQ.1) -AND.(J3.EQ.1)) GO TO 16 

IF (MP1.EQ.1). GO TO 16 

IF(NP1.EQ.1). GO TO 16 

IF (IFBP1.GT.NBP1) GO TO 16 

IF(MP1.£Q.IRP1) GO TO 44 

IF(NP1.EQ.IRP1) GO TO 44 

IF(IFAP1.NE.1) GO TO 44 

IF(IFBP1,EQ.IFAP1) GO TO 16 

IF(IFBP1.GT.NBP1) GO TO 16 4, © 
IFAPFB=LFAPIt+tIFBP1-2. a 
IF(IFAPFB.LE.1) GO TO 24 
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RULE = 8./IFAPFB 
DIFPRO= ((IFBP1-1.)/ (NBP1-1) )-((IFAP1-1.) /(NAP 1-1) ) 
DIFNEG=-DIFPRO 
24 £F(03.E0.1}) GO’ To 23 
IF(NP1.EQ.IRP1) GO TO 26 
IF (IFBP1.EQ.1) GO TO 23 
IF(IFAPFB.LE.1) GO TO 26 
IF(DIFPRO.LT. RULE) GO TO 31 
R (NP 1, IFAP 1,33) =0. 


GO TO 23 
an IF(DIFNEG.GT.RULE) GO TO 23 
26 R (NP1,I1PAP1,33) =P*R1(NP1, IFAP1,J3) +Q*S1(NP1,IFAP1+1,J3-1) 


IF(R(NP1,IFAP1,J3) .LT.10E-10) R(NP1,I1FAP1,J3) =0. 
23 IF(J3.EQ.3) GO TO 16 

IF (MP1.EQ.IRP1) GO TO 27 

IF(IFAP1.EQ.1) GO TO 16 

IF (IFAPFB.LE.1) GO TO 27 

IPF(DIFPRO.LT.RULE) GO TO 33 

S (NP1,IFAP1,J3) =0. 


GO TO 16 
33 IF (DIFNEG.GT. RULE) GO TO 16 
a S (NP 1, [FAP 1,33) =P1*S1(NP1-1,IFAP1,J3) + 
a Q1*R1(NP1-1,IPAP1,J3+1) + 1. 
IF(S(NP1,IFAP1,J3) .LT.10E-10) S(NP1,IPAP1,J 3) =0. 
16 R1(NP1,IPAP 1,J3) =R (NP1,IFAP1,J3) 
ss S$ 1 (NP1,IFAP1,33)=S (NP1, IFAP1,J3) 
15. CONTINUE 
AC ENB GIVES THE EXPECTED NUMBER OF OBSERVATIONS ON THE 
Cc POORER TREATMENT. 


ENB= (R1(IRP1, 1,2) +S1(IRP1,1,2)) /2- 
P3=(P+P1)/2. 
WRITE (5,210) P3,ENB 
210 FORMAT('0',T52,F8.4,T72,F12.6) 
STOP 
END 
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